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This article reviews the results of a test of 
the Gibbs-DiMarzio theory by Monte 
Carlo Simulation. The simulation employed 
the bond-fluctuation model on a simple 
cubic lattice. This model incorporates two 
kinds of interactions: the excluded 
volume interaction among all monomers of 
the melt and an internal energy of the 
chains, which favors large bonds and makes 
the chains stiffen with decreasing 
temperature. The stiffening of the chains 
leads to an increase of their volume 
requirements, which competes with the 
packing constraints at low temperatures. 
This competition strongly slows down the 
structural relaxation of the melt and 
induces the glassy behavior. The model 
therefore takes into account the main 
opposing forces which the Gibbs-DiMarzio 
theory makes responsible for the glass 
transition. For this model the entropy was 



calculated from the internal and the free 
energy (derived from the chemical potential 
and the single chain partition function) 
and compared with various theoretical 
predictions: the Gibbs-DiMarzio theory, a 
theory by Flory for semiflexible polymers 
and an extended theory by Wittmann 
considering Milchev's criticism on Flory's 
calculation. The latter extended theory 
provides the best description of the 
simulation data. 
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bond vector of length h 
specific heat 
constant 

diffusion coefficient at infinite tempera- 
ture 

diffusion coefficient at temperature T 
energy density as a function of monomer 
density p and reciprocal temperature j5 
energy difference in the Metropolis crite- 
rion 

free energy density as a function of 
monomer density p and reciprocal temper- 
ature j5 



f flexibility parameter [see Eq. (9)] 

H number of holes (empty lattice sites or unit 

cells) 

df€{b ) bond vector energy function of the bond- 

fluctuation model 

^B Boltzmann constant 

K number of polymers per simulation cell 

L linear dimension of the simulation cell 

M total number of lattice sites 

N chain length 

A^o number of overlaps between a test chain 

and the melt [see Eq. (7)] 
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approximate number of adjacent empty 
lattice sites for placing of a new monomer 
during the construction of a polymer [see 
Eq. (12)] 

Flory's approximation for A^empty [see Eq. 
(13)] 

Gibbs-DiMarzio's approximation for 
A^empty [scc Eq. (14)] 

Milchev's approximation for A^empty [see 
Eq. (15)] 

addition probability for a bond vector b 
[see Eq. (4)] 

insertion probability for a test chain in a 
melt at monomer density p and reciprocal 
temperature jB [see Eq. (6)] 
wave number 

entropy density as a function of monomer 
density p and reciprocal temperature jB 
Flory's approximation for the entropy den- 
sity [see Eq. (16)] 

Gibbs-DiMarzio's approximation for the 
entropy [see Eq. (17)] 
Milchev's approximation for the entropy 
[see Eqs. (16) and (17)] 
collective static structure factor 
temperature 

critical temperature of mode-coupling 
theory 

(kinetic) glass transition temperature 
Kauzmann temperature 
Vogel-Fulcher temperature 
potential experienced by a bond when 
added to a chain of (N-l) monomers [see 
Eq. (4)] 

potential between an inserted test chain (t) 
and the melt (m) 

soft potential between an inserted test 
chain (t) and the melt (m) [see Eq. (7)] 
coordination number of the lattice 
single chain partition function [see Eq. (5)] 

\/T 

energy parameter of ^(b) (units: e/ 

^B=l) 

energy parameter of the bond angle in the 

Gibbs-DiMarzio and the other theories 

[see Eq. (9)] 

effective energy barrier (in temperature 

units) of the relaxation time [see Eq. (22)] 

isothermal compressibility 

excluded volume parameter [see Eq. (7)] 

excess chemical potential [see Eq. (6)] 



Vk number of configurations for the /:th chain 

if there are already (k-l) chains on the 

lattice [see Eq. (12)] 
p monomer density (= KN/M) 

pGff(T) effective monomer density (= K N/ 

[KN-\-H]) [seeEq. (20)] 
T Structural relaxation time 

cf) volume fraction of occupied lattice site in 

the bond-fluctuation model (=SKN/L^) 
fl total number of melt configurations [see 

Eq. (8)] 
i^intra intrachain contribution to fl [see Eq. (10)] 

i7inter intcrchaln contribution to fl [see Eq. (1 1)] 

1, Introduction 

In 1956 Julian H. Gibbs wrote a seminal paper [1], in 
which he suggested that a recent theory by Flory [2] 
could be applied to explain the glass transition of 
polymer melts. In his study Flory obtained a negative 
configurational entropy for densely packed polymers of 
sufficient length if their stiffness exceeded a critical 
value at low temperatures. Flory avoided this paradox 
by invoking crystallization, whereas Gibbs pointed out 
that most polymers due to their structural complexity 
rather undergo a glass transition. Therefore he inter- 
preted Flory' s result as a possible explanation of 
the Kauzmann paradox [3-5], i.e., of the experimental 
observation that the configurational entropy difference 
between the (supercooled) liquid and the crystalline 
state seems to vanish at a finite temperature in the 
limit of an infinitely slow cooHng process. An extensive 
recent survey of experimental data [6,7] suggests that 
the extrapolated temperature of the vanishing entropy 
difference, the Kauzmann temperature Tk (< Tg), 
roughly coincides with the Vogel-Fulcher temperature 
Tq resulting from fitting the Vogel-Fulcher equation 
[see Eq. (1)] to pertinent relaxation times or to the 
viscosity. 

This interpretation of Gibbs motivated the devel- 
opment of the Gibbs-DiMarzio theory [8,9], which 
views the experimental glass transition as a signature of 
an underlying thermodynamic (second-order) phase 
transition. This phase transition occurs at Tk and is 
driven by the conflict between the bulky low-energy 
polymer conformations (increasing stiffness) and the 
contraction of the sample's volume at low temperatures. 
At high temperatures the polymers should be able to 
easily reconcile these opposing forces so that the 
entropy per monomer s is larger than zero. When the 
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temperature is reduced the mentioned conflict gradually 
starts to develop and decreases the entropy until the 
number of accessible configurations has become so 
small that s vanishes in the thermodynamic limit. The 
corresponding state of the melt is called the "ground 
state of amorphous packing," in which the melt remains 
for T<Tk. 

From this general theoretical treatment many specific 
predictions have been derived and tested in experiments, 
such as the variation of the glass transition temperature 
Tg (assuming that T^ behaves as T^) with pressure [10] 
or with molecular weight [8,11-14], the temperature 
dependence of the specific heat at Tg [15,16], the influ- 
ence of crosslinks [17-19] and plasticizer [20] on the 
glassy behavior of the melt, the composition dependence 
of Tg for copolymers and polymer mixtures [21], etc. 
These studies show that the Gibbs-DiMarzio theory 
provides a rather good description of the experimental 
situation. 

However, to the best of our knowledge, in all tests of 
the theory only derived quantities were used. No attempt 
has been made to compare the temperature dependence 
of the basic quantity, the configurational entropy (which 
is admittedly hard to measure), with experimental data. 
With the present work we want to present such test for 
a computer model of a glassy polymer melt. 

This model and some of its properties are discussed 
in the Sec. 2. Section 3 then gives an outline of the 
method used to determine the entropy, whereas the 
Sec. 4 exposes those aspects of the Gibbs-DiMarzio and 
of related theories, which are necessary for the sub- 
sequent comparison. This comparison is discussed in 
Sec. 5, and finally Sec. 6 contains our conclusions. 



2. Definition and Properties of the Model 
2.1 The Bond-Fluctuation Model for Glassy 
Polymer Melts 

The bond-fluctuation model [22] was proposed as an 
intermediate description between a highly flexible 
continuum treatment and the traditional lattice models, 
where the flexibility of the chains, described as random 
or self-avoiding walks [23], is fully determined by the 
structure of the underlying lattice. It shares with the 
latter models the simple — and from a computational 
point of view highly efficient [24] — lattice structure, but 
distinguishes itself from them by exhibiting a multitude 
of bond vectors, as also present in continuum models. 
Due to the lattice structure a monomer of the model 
does not directly correspond to a chemical monomer, 
but rather to a group of chemical monomers (typically, 
such a group comprises 3-5 monomers for simple 
polymers, such as polyethylene; see Fig. 1 [25]). Since 
a lattice bond should thus be interpreted as the vector 
joining the mass centers of these groups, its length and 
direction will fluctuate. The bond-fluctuation model 
maps this idea onto a square or simple cubic lattice by 
associating a monomer with a unit cell of the respective 
lattice. In order to impose local self-avoidance of the 
monomers (excluded volume interaction) and to guaran- 
tee uncross ability of the bond vectors during the course 
of the simulation (no phantom chains) the allowed bond 
vectors are obtained from the sets {(2,0,0), (2,1,0), 
(2,1,1), (2,2,1), (3,0,0), (3,1,0)} by all symmetry 
operations of the lattice. This yields 1 08 different bond 
vectors. 




^ 




Fig. 1. Schematic illustration of the construction of a coarse-grained model for a macromolecule such as 
polyethylene. In the example shown here, the subchain formed by the three C-C bonds labeled 1,2,3 is 
represented by the effective bond labeled as I, the subchain formed by the three bonds 4,5,6 is represented by 
the effective bond labeled as II, etc. In the bond-fluctuation model the length b of the effective bond is allowed 
to fluctuate in a certain range fomin ^ b ^ b^^, and excluded-volume interactions are modeled by assuming that 
each bond occupies a plaquette (or cube) of 4 (8) neighboring lattice sites which then are all blocked for further 
occupation. From Ref. [27]. 
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In addition to the excluded volume interaction an 
energy function W(b) is associated with each bond 
vector b , which favors bonds of length b = 3 and direc- 
tions along the lattice axes (i.e., ^(b) = 0) in compari- 
son to all other available bond vectors (i.e., ^(b) = t) 
[26-28]. Figure 2 illustrates the effect of this energy 
function. When the temperature decreases each bond 
tries to reach the ground state (i.e., a bond with 
^(b) = 0) and thereby blocks four lattice sites for other 
monomers. This loss of available volume generates a 
competition between the energetically driven expansion 
of a bond and the packing constraints of the melt. Due 
to this competition some bonds are forced to remain in 
the excited state (see Fig. 2). They are geometrically 
frustrated [26-28]. The development of the geometric 
frustration during the cooling process causes the glassy 
behavior of the model and corresponds to the driving 
force which the Gibbs-DiMarzio theory makes respon- 
sible for the glass transition of polymer melts [8]. 

For the present study a chain length of A/^= 10 was 
used. If one takes into account that a lattice monomer 
roughly corresponds to a group of three to five chemical 
monomers (see above) our simulation deals with fairly 
short, oligomeric chains. The cubic simulation box (of 
linear dimension L = 30) contained K= \S0 chains so 
that the volume fraction of occupied lattice sites is 



cl)=SNK/L^= 0.53. This value is a compromise 
between two requirements: it is high enough for the 
model to exhibit the typical behavior of dense melts at 
high temperatures [24] as well as pronounced frustration 
effects at low temperatures [26-28] and low enough to 
allow for a sufficient acceptance rate of monomer (or 
chain) moves to make the equilibration of the melt in the 
interesting temperature range possible [29]. In order to 
improve the statistics 16 independent simulation boxes 
were treated in parallel. Thus the total statistical effort 
involves 28800 monomers, which ensures a high statisti- 
cal accuracy of the results. 

Due to the above mentioned development of the 
geometric frustration with decreasing temperature the 
structural relaxation time strongly increases if the usual 
bond-fluctuation dynamics is used. This makes the 
equilibration of the melt very time consuming at low 
temperatures. The usual bond-fluctuation dynamics 
consists in moving a randomly chosen monomer in a 
random direction along the lattice axes. These local 
moves are supposed to mimic a random force exerted on 
a monomer by its environment. They lead to Rouse-like 
dynamics which is typical of a polymer in a dense melt 
[30] . However, since we are only interested in 
equilibrium properties of the melt, this realistic dynam- 
ics may be replaced by an artificial one which 




= 6 



- — I — 1 — ' — — '*^ — I — w- 



jump forbidden 



Fig. 2. Sketch of a possible configuration of monomers belonging to different chains in the 
melt in order to illustrate the effect of the model's energy function and the concept of 
geometric frustration. All bond vectors shown in this picture have the energy € except the 
vector (3,0,0) which is in the ground state. This vector blocks four lattice sites (marked by o) 
which are no longer available to other monomers, since two monomers may not overlap. Due 
to the excluded volume interaction the jump in direction of the large arrow is also forbidden. 
This leads to geometric frustration. From Ref. [27]. 
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uses nonlocal moves to speed up the equihbration. 
A nonlocal move involves a collective motion of all 
monomers of a chain. Such a collective motion may be 
realized by the so-called slithering-snake dynamics 
[23,31]. for instance. In the slithering-snake dynamics 
one attempts to attach a randomly chosen bond vector 
(from the set of allowed bond vectors) to one of the ends 
of a polymer (both also randomly chosen). If the attempt 
does not violate the excluded volume restriction, the 
move is accepted with probability exp [-AE/kBT], 
where AE" is the energy difference between the newly 
added bond and the last bond of the other end of the 
chain. This last bond is removed, provided the move is 
accepted. By means of the probability exp [-AE/kBT] 
(Metropolis criterion [32]) temperature is introduced 
into the simulation [33]. 

2.2 A Brief Review of the Model's Properties 

The above mentioned competition between the inter- 
nal energy and the packing constraints of the polymers 
prevents the melt from crystallizing so that it may easily 
be supercooled. During the supercooling process the 
chains stiffen [26,27,29] and the structural relaxation 
time of the melt increases [26,27]. 



This increase of the structural relaxation time is 
exemplified in Fig. 3 which shows the temperature de- 
pendence of the polymer diffusion coefficient [27]. In 
the available temperature interval the diffusion coeffi- 
cient decreases by about two orders of magnitude and its 
temperature dependence may be rather well fitted by a 
Vogel-Fulcher equation 



D(T) = D^cxp 



T-To 



(1) 



yielding /)«> = 8.61 X 10"" ± 0.32 X 10-^ C = 0.396 ± 
0.041 and To = 0.17 ± 0.02. However, the resulting 
Vogel-Fulcher temperature To — 0.17 certainly over- 
estimates the absolute freezing temperature of the 
model considerably, since the extrapolation is based on 
a variation of the diffusion coefficient of two decades 
only, and one knows from experiments [34-36] where 
ten decades are available that the fitted values of this 
Vogel-Fulcher temperature depends on the dynamic 
range available. The extrapolated result for the Vogel- 
Fulcher temperature should therefore not be looked 
upon as an accurate value for the absolute freezing point 
of the studied model, but rather as an estimate for the 
crossover temperature to the regime in which the 
interesting glass physics of this model occurs. 



0.0001 : 



le-05 : 




Fig. 3. Plot of the diffusion coefficient versus temperature.The dashed line corresponds to a fit by the Vogel- 
Fulcher equation [Eq. (1)]. From Ref. [27]. 
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Due to the strong slowing down of the mobility of 
the polymers the melt starts to vitrify as soon as its 
structural relaxation time matches the time scale of the 
simulation. The resulting (kinetic) glass transition 
temperature depends on both cooling rate [26,27] and 
chain length [26,37]. It varies with the inverse chain 
length in a linear and with the cooling rate in a non- 
linear fashion. Both of these results are also found in 
experiments [4,38]. 

During the vitrification process the melt maintains its 
amorphous structure. This is illustrated in Fig. 4 which 
depicts the wave number (q-) dependence of the collec- 
tive structure factor of the melt at two representative 
temperatures, one characteristic of the equilibrium 
liquid state (T= 1.8) and one characteristic of the frozen 
glassy state (r=0.05) [26,28]. The structure factor 
exhibits the shape expected for an amorphous material. 
It is small at small q -values (where it coincides with the 
isothermal compressibility), then increases and devel- 
ops a maximum (the so-called "amorphous halo" [39]) 
at the typical distance between neighboring monomers 
(2-3 lattice constants) before it decreases again. This 
behavior remains unchanged with decreasing tempera- 
ture. The strong temperature dependence for q>3J 
should not be interpreted physically, since these large 



^-values probe the length scale of a lattice constant, 
where artifacts due to the underlying lattice must 
become visible. The important feature of the structure 
factor is the variation at small and intermediate 
^-values, where it reproduces well-known experimental 
results [39]. 

If one further equilibrates the melt in the temperature 
region close to the kinetic glass transition the model 
exhibits the phenomenon of physical aging [4,40]. A 
first analysis of these effects for the present model [41] 
shows that the approach towards equilibrium is an 
extremely stretched process which obeys an aging-time- 
temperature superposition property. Similar results are 
also found in experiments [4,40]. 

If one removes these nonequilibrium effects 
sufficiently the incoherent intermediate scattering 
function decays in two steps [42,43]. This two-step 
relaxation behavior is predicted by the mode-coupling 
approach [43-45] to the structural glass transition and 
can be described quantitatively by the theory. For the 
critical temperature of the theory this analysis yields a 
value of Tc ~ 0.15. Since experimentally Tc > To [43], 
this result is another evidence that the above shown 
Vogel-Fulcher fit overestimates the absolute freezing 
temperature of the model. This conclusion is also 





1.6 


+ 




1.4 


T 


sym. 




0.05 


** 






Kt/K? 0.05 


+ + 




1.2 


1.80 


^ o ^ 






Kj/K° 1.80 


♦ +++0 . ++ 


_, 


1 


r % **-*:-^^^ *\ - 


u- 




^* 




0.8 


o 


CO 




\ 




0.6 






0.4 






0.2 


.^^^+.*.+.t.ttt «. ^ , 




1 1 1 1 1 







Fig. 4. Plot of the collective static structure factor 5 (^) vs ^for r= 0.05 (O; glassy phase) and T= 1.8 (+; liquid 
phase). The cooling rate used in the simulation was Fq = 4 X 10~^. Additionally, a dashed and a dotted horizontal 
line are shown which are the isothermal compressibilities kt/kt= TpKr at the respective temperatures. From 
Ref. [28]. 
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corroborated by a recent extensive analysis of many 
different dynamic properties (relaxation function of the 
bond vectors, of the radius of gyration, etc., mean- 
square displacements of the monomers and chains, 
Rouse-mode analysis), which were calculated from 
completely equilibrated configurations in an extended 
temperature regime close to the previously determined 
To [46]. This study yields To ^ 0.12-0.13 < T^, as in 
experiments. 

This keyword-like summary of important properties 
of the model should illustrate that despite its simplicity 
the model reproduces experimental findings of (fragile 
[6]) structural glass formers. Therefore it is a good 
working model to determine the entropy and to compare 
the results with the predictions of the Gibbs-DiMarzio 
theory. 



3. Entropy and Monte Carlo Simulation 

This section briefly describes the main steps and 
quantities, which are used to determine the entropy by 
Monte Carlo simulation. To this end, consider a (simple 
cubic) lattice with M sites and K polymers. Then the 
monomer density p is given by p = N K/M. In order to 
calculate the entropy density s(p, jS) one has to know 
the energy and the free energy density as a function of 
p and reciprocal temperature jS , since 



s(p,l3) = l3e(p,l3)-mp,P) 



(2) 



Whereas the energy density can easily be measured in 
the simulation as an average of the model's energy 
function (see Sec. 2), the free energy density is more 
difficult to obtain. It may be expressed as [47,51] 



WCp.W'^ 



,„ £-■-,„ *'*'*'"■ 



A^ 



A^ 



-r 

NJo 



dp>exfp'jS), 



(3) 



where ^p is the partition function of a single polymer 
and ^Ltex is the excess chemical potential. 

In order to calculate ^p we used a method proposed 
by Kumar et al. [48], whose basic idea is as follows: 
Imagine a chain of (A^-1) monomers, to which an 
additional bond vector shall be attached. This bond 
vector Z» AT.] experiences a potential V(bn..i) which deter- 
mines the addition probability /7(Z?Ar_i) by 

^ ' , o. = E <exp [-jsy(Z^^_i)]) = m%p{b^.,) . 

=^pVV l,p>l 6jy_i=i 

(4) 



Iteration of this equation yields 

3£^{N, P) = M 108"-' p{b,)...p (6„_,) . 



(5) 



Therefore the desired partition function may be ex- 
pressed in terms of addition probabilities for successive 
bonds which are accessible in the simulation [47,49,50]. 
The second input quantity for the calculation of the 
free energy is the excess chemical potential. Physically 
the chemical potential is the change in free energy if a 
test chain (i.e., an additional chain) is added to a 
polymer melt. When inserted in the melt (m) the test 
chain (t) experiences a potential Vmt with the surround- 
ing polymers, the Boltzmann factor of which, averaged 
over all melt and test chain configurations, determines 
the insertion probability p-^^^ and the excess chemical 
potential by [51] 



I3p...(p,l3) = - In <exp [-jS V^t]) = - In pUp.P) ■ 



(6) 



Quite generally, the insertion probability for polymers 
is a very small quantity. For instance, for the present 
model it is /7ins ~ 10"^^ already at infinite temperature 
and decreases still further with falling temperature [47]. 
The reason for these small values is that a randomly 
chosen test chain configuration, inserted at random in 
the melt, is likely to overlap many times with the 
surrounding polymers. Since every single overlap 
forbids the insertion, successful attempts stem from the 
events that the chosen configuration happens to fit in 
the cavity of the melt at the point of insertion. Since 
these events are very seldom, the random insertion 
method, which can successfully be applied to small 
molecules, is very inefficient for polymers. In order to 
circumvent this inefficiency various methods have been 
designed [49-51]. We use a method proposed by Miiller 
et al. [49,50], whose basic idea is as follows: Instead of 
strictly forbidding overlaps between the test chain and 
the melt every overlap is penalized by a finite potential 



-l3Vrm(No,X)=No(x^^)\nX, 



(7) 



where A^o(^mt) is the number of overlaps, which depends 
on the melt and the test chain configuration x^t, and A 
is a tunable parameter ranging between and 1 . The 
value A = corresponds to complete excluded volume 
interaction, whereas A = 1 means no excluded volume 
interaction. The intermediate A -values between these 
two extreme states are used as auxiliary variables in the 
(extended-ensemble) simulation in order to gradually 
approach the desired case of complete excluded volume 
interaction and to gather sufficient statistics for it. 
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4. Entropy Approximations: 

Gibbs-DiMarzio and Related Theories 

In the following section the simulation results for the 
entropy will be compared with the prediction of the 
Gibbs-DiMarzio theory [8], of a theory by Flory [2] and 
of a theory by Wittmann [52] which considers 
Milchev's criticism on Rory's calculation [53]. The 
present section is devoted to introduce those aspects of 
these theories which are needed for the comparison. 

Consider a lattice with M sites and K polymers of 
length A^. Each monomer is supposed to occupy a single 
lattice site (contrary to the bond-fluctuation model; see 
Sec. 2). Furthermore, let £2 denote the total number of 
configurations of the melt. All of the above mentioned 
theoretical approaches factorize £2 into an intra- and an 
interchain contribution so that the entropy density is 
approximated by 



:lln/2«lln[a 



a^ 'ir 



(8) 



The intrachain contribution accounts for the increase 
of the chain's stiffness during supercooling. This is 
modeled in the theories by associating an energy with 
the bond angle. The energy is if two successive bonds 
are colinear and g otherwise. For such a two-level 
system the probability that two successive bonds are not 
cohnear ("flexed"), the so-cMcd flexibility f, reads 



/= 



(z-2)exp[-;8g] 
l+(z-2)exp[-|8e] 



(9) 



where z is the coordination number of tiie lattice. With 
this assumption about the internal energy of a polymer 
the intrachain contribution to O is given by the follow- 
ing binomial distribution [52] 



n. 



\fK(N-2)/ \z-l 



^(K(N-2)\(J_Y-"^'"-'' (z:^ 



"■ \) .(10) 



In addition to ilintra an expression is needed for the 
number of different ways in which K polymers may be 
placed onto the lattice. This interchain contribution to Q 
is given by 



^-1 

a - — — n 

i Winter — ^K t^ y 11 



Vk+l 



(11) 



It involves the total number of configurations Vk+i that 
the (^-1-1 )th chain may adopt if there are already k chains 
on the lattice. Approximately, Vk^i consists of three 
factors 



V,,, - (M-W)XiV^-;,Xz(z- 1)^ 



(12) 



The first factor is the number of remaining empty lattice 
sites and thus of potential starting points for the (^-i-l)th 
polymer, whereas the last factor represents the number 
of possibilities to place the remaining (A^- 1) monomers 
onto the lattice if immediate backfolding is forbidden. 
The second factor approximately accounts for the 
self- and mutual avoidance of the polymers. It is related 
to the probability that the newly chosen site for the next 
monomer of chain (^+1) is in fact empty. The three 
theories mentioned above differ from each other in the 
ansatz made for A^empty Flory uses the volume fraction of 
empty lattice sites for A^empty, ie.. 



M 



empty, F " 



kN 
~M ' 



(13) 



whereas Gibbs and DiMarzio and Milchev work with 



N. 



1-kN/M 



empty, GDM - j _^ ^^_ l)/(^^2) 



and 



N, 



l-kN/M M-kN 



empty, M j_^/^ M-kM/K ' 



(14) 



(15) 



respectively. They divide Flory's ansatz by expressions 
which approximately take into account that not all 
empty lattice sites can serve as starting points for the 
new polymer, but only those which lie outside of the 
volume already consumed by the other k chains. To 
estimate this volume Gibbs and DiMarzio argue that 
there are Mz/2 different nearest neighbors (on a 
lattice with periodic boundary conditions) so that 
(N-l)/(Mz/2) is the volume fraction occupied by a 
chain, whereas Milchev assumes that a chain consumes 
on average M/K sites. 

With these different ansatzes the entropy density can 
be calculated. The final results read [52] 



^f(/ K, HH,z) = suif. K, H H, z)- 



in Flory's approximation and 



KN 
KN+H 



[-;] 



(16) 



ioDM(/ K,N,H,z) = 

rfZr^lTJ ^ zjKN+H) , (KN+H) 



(17) 



for the Gibbs-DiMarzio theory. In Eqs. (16) and (17) ^m 
is the results of Milchev's theory, H denotes the number 
of holes (i.e., of empty lattice sites: H = M-KN) and 
N is number of nearest neighbors of a (rodlike) 
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polymer [52]. In the limit of the completely filled lattice 
{H -> 0), of A^ -> 00 and r -> (i.e.,/ -> 0), Milchev's 
entropy vanishes [52] so that 



H,f^O,N- 



and 



^f(/ K K H, z) 



SGDuif, K, N, H, z) > ^ 



(18) 



X 



r 2i 




r 2i 


I-- 


In 


I-- 


z_ 




z_ 



->-l 



(19) 



The negative entropy of Eq. (19) in the low temperature 
limit is interpreted as the theoretical counterpart of the 
Kauzmann paradox [3-5]. The Gibbs-DiMarzio theory 
resolves this paradox by associating a second-order 
phase transition from an ideal glass to a liquid phase 
with the temperature Tk at which the configurational 
entropy vanishes. Such a transition does not occur in 
Milchev's theory, where the entropy remains positive 
for any finite temperature [52]. This means that the 
choices made in Eqs. (13)-(15) do not only introduce 
quantitative corrections to the final result, but lead to a 
qualitatively different behavior. It is therefore interesting 
to test which of the theoretical approaches describes the 
simulation data better. 

5. Comparison of Theory and Simulation 
5.1 Determination of Input Parameters 

In order to perform a comparison between the above 
outlined theories and the simulation data the relevant 
input parameters, the flexibility/, the number of holes 
H and the coordination number z of the lattice, have to 
be determined for the bond-fluctuation model. 

The theories associate a two-level system with the 
bond angle, whereas the present model uses a two-level 
system for the bond length. Therefore we replace 
the flexibility by the fraction of bonds in the excited 
state. The temperature dependence of the so-defined 
/-parameter is shown in Fig. 5 [54]. At high tempera- 
tures almost all bonds are in the excited state. As the 
temperature decreases more and more bonds enter the 
expanded ground state (see Fig. 2) and the /-parameter 
decreases. In addition, the expansion on the length scale 
of a bond also induces an expansion of the chains as a 
whole, which entails an increase of the chain's stiffness 
[29]. This is exactly the effect that the theories try to 
model by the introduction of the flexibility. Note that 
simulations using two-level systems associated with the 
bond angle (as in the analytic theories) give evidence 
for the formation of strong nematic-type short-range 
order among the chains, which is not typical of the glass 
transition of most polymers, and hence such models are 



less suitable for testing the entropy theory of the glass 
transition. 
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Fig. 5. Fraction/of bonds in the excited state versus inverse temper- 
ature. FromRef. [54]. 

The second input quantity is the number of holes. In 
order to determine H remember that a monomer of the 
bond-fluctuation model does not correspond to a single 
lattice point (as in the theories), but to a whole unit cell 
of the simple cubic lattice. Therefore a hole should be 
interpreted as an empty unit cell. With this identification 
one can define an effective density peff by [54] 



Peff = 



^A^ 



KN + H{T)' 



(20) 



The temperature dependence of this density is shown in 
Fig. 6 [54]. For all temperatures peff is larger than the 
volume fraction of occupied lattice sites (f) = 0.53, which 
illustrates that the melt is in fact much denser than could 
have been anticipated on the basis of cf) alone. Further- 
more, the effective density increases with decreasing 
temperature, since the accompanying increase of the 
number of bonds in the ground state reduces the avail- 
able volume and thus makes the density rise. 

Finally, the coordination number z of the lattice has to 
be determined. In the theories, z coincides with the 
number of nearest neighbors. Therefore a natural choice 
is to take z as the number of monomers in a sphere with 
radius b^^ = Vio (see Sec. 2) around a central 
monomer. Figure 7 shows that the resulting z -values 
range between 11 and 12, which is typical of dense 
(amorphous) packing and that z decreases with decreas- 
ing temperature due to the expansion of the bond vectors 
[54]. 
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Fig. 6. Effective density peff versus inverse temperature. From Ref. 
[54]. 
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Fig. 8. Comparison of the temperature dependence of the entropy 
per lattice site s with the Milchev, Flory, and Gibbs-DiMarzio theory 
[see Eqs. (16)-(17). From Ref. [54]. 
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Fig. 7. Coordination number z versus temperature. From Ref. [54]. 



5.2 Results 

Having made the above discussed identifications,/, H 
and z can be inserted in Eqs. (16)-(17). The resulting 
temperature dependence is compared with that of the 
simulation data in Fig. 8 [54]. As expected, the 
simulated entropy decreases with decreasing tempera- 
ture. However, instead of exhibiting a clear tendency to 
vanish at a finite temperature its decrease gradually 
weakens with falling temperature so that it rather seems 



to stay larger than zero also at very low temperatures, 
where no further simulation points are up to now 
available. Therefore IMilchev's theory lies closer to the 
simulation data than the Gibbs-DilMarzio or the Flory 
theory. 

However, taking into account that Tk = in Milchev's 
theory [52] and that there is the Adam-Gibbs relation 
[55] 



D{T) = D^ exp 



I Ts(T)\' 



(21) 



which was proposed as a possible rationalization of the 
Vogel-Fulcher equation, the found temperature depen- 
dence of the entropy challenges the Vogel-Fulcher 
analysis of the diffusion data (see Fig. 3) [54]. 
Therefore Fig. 9 directly compares Eq. (21) with the 
measured diffusion coefficients of the chains. Since 
there are only two data points for the entropy at r= 0.3 
and r=0.6 in the temperature interval where D is 
available, the prefactor in Eq. (21) was taken from the 
simulation (Da= = 5.2X 10"^) and the constant C was 
adjusted in such a way that Eq. (21) either coincided at 
T= 0.3 or at r= 0.6 with the simulation data. Despite 
the preliminary status of this comparison the result for 
r=0.3 is not unreasonable. This suggests that the 
Adam-Gibbs formula could be a valid description for 
the temperature dependence of the structural relaxation 
times of the present model and that the absolute 
freezing point should be smaller than To — 0.17. 
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Fig. 9. Comparison of the temperature dependence of the diffusion 
coefficient with the Adam-Gibbs formula [Eq. (21)]. See text for 
further details. From Ref. [54]. 



6. Concluding Remarks 

This paper summarizes the results about the tempera- 
ture behavior of the entropy for a simple model of a 
supercooled polymer melt. The model consists of linear 
monodisperse chains on a (simple cubic) lattice, which 
interact by a hard core potential and possess an internal 
energy that increases the chains' stiffness with decreas- 
ing temperature. During supercooling a competition 
therefore arises between the tendency of a single chain 
to expand and the corresponding enlarged volume 
requirement of all chains. These opposing forces 
strongly slow down the structural relaxation and mimic 
the mechanism that the Gibbs-DilVIarzio theory suggests 
to be responsible for the glassy behavior of polymer 
melts. In fact, this suggestion partly inspired the design 
of the model. 

However, when the temperature dependence of the 
entropy is compared with the theoretical prediction one 
finds that the shape of the curves are very similar, but 
that the simulation data are shifted upwards with respect 
to the theory. A similar shift is also obtained theoreti- 
cally by Wittmann's analysis of IVIilchev's criticism. 
Therefore IVIilchev's theory approximates the simulation 
data better than Gibbs-DilVIarzio's or Flory's theory. 
This finding leads to several conclusions and questions: 
(1) Presumably, the Gibbs-DilVIarzio theory (as well as 
the other theories) can well describe the temperature 
dependence of the model's specific heat c, since 



c = Tds/dT, as it is also found in experiments. 
In order to test this conjecture more data points are 
needed at low temperatures. This is a prospective topic 
for future work. (2) Gibbs-DilVIarzio's correction of 
Flory's ansatz for A^empty [see Eq. (14)] only implies 
minor changes for the entropy. Both theories predict the 
entropy to vanish at about the same temperature (i.e., 
Tk ~ 0.17 -0.18). Interestingly, this temperature coin- 
cides with the Vogel-Fulcher temperature To obtained 
from the diffusion coefficient. Whether this is a mere 
coincidence or whether there is an explanation on the 
basis of the model's properties is still unclear. (3) The 
steepness of the entropy decrease in the simulation 
becomes weaker with decreasing temperature. This 
tendency to rather level off than to vanish suggests that 
the entropy remains larger than zero also at very low 
temperatures, where no data points are up to now 
available. Certainly, further simulations are necessary to 
clarify this behavior. (4) An interesting question in this 
context is the validity of the Adam-Gibbs relation [see 
Eq. (21)]. This relation implies for the relaxation 
time T ^ cxp[A/Ts], contrary to the usual expectation 
T ^ cxpljBf] ^ exp[^/A:B]. Our simulation only 
provides a first indication that the Adam-Gibbs relation 
could be applied. To test this indication the structural 
relaxation time has to be determined at lower tempera- 
tures. Work in this direction is underway [46]. Since 
IVIilchev's theory is closest to the simulation data, it 
would then also be possible to check its prediction for 
the temperature dependence of the relaxation time, i.e., 
[56] 



T(r)ocexp 



[(1)1 



2c 
with a = — 

z 



(22) 



In order to apply this equation an extension of entropy 
measurement is also needed due to the dependence of 
the exponent a on the specific heat. 
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